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Abstract. A construction of "sparse potentials", suggested in 
[10] for the lattice 1 d : d > 2, is extended to a wide class of com- 
binatorial and metric graphs whose global dimension is a number 
D > 2. For the Schrodinger operator —A — aV on such graphs, 
with a sparse potential V, we study the behavior (as a — > oo) of 
the number A_(— A — aV) of negative eigenvalues of —A — aV. 
We show that by means of sparse potentials one can realize any 
prescribed asymptotic behavior of A_(— A — aV) under very mild 
regularity assumptions. A similar construction works also for the 
lattice Z 2 , where D = 2. 



1. Introduction 

There exists a far-reaching parallelism between the theory of the 
Schrodinger operator on M d and its discrete analogue on Z d . Still, there 
are issues where this parallelism is violated. Certainly, the most known 
examples of such violation come from the fact that, unlike the case of 
M. d , the Laplacian on Z d is a bounded operator. However, there are 
also examples of a different origin. One of them concerns estimates on 
the number AL(— A — aV) of the negative eigenvalues for the operator 
—A — aV, where V > is a potential and a > is a large parameter 
("the coupling constant"). Both for M. d and for Z d the "Rozenblum - 
Lieb - Cwikel estimate" (RLC estimate) is satisfied: namely, one has 

(1.1) AL(-A - aV) < Ca^WVfJ*, C = (7(d); d > 3. 

Here || ■ \\ q stands for the L 9 -norm in the continuous case and for the 
£ 9 -norm j n the discrete case. In both cases the inequality (1.1) is valid 
for any potential V such that the norm in the right-hand side is finite. 
The difference between these cases is this: 

on R d : i\L(-A - aV) = o(a d/2 ) if and only if V = 0, 
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while 

(1.2) onZ d : AL(-A - aV) = o{a d/2 ) for any V G £ d ^ 2 . 

This latter effect was, probably, observed for the first time in [9]. The 
same effect manifests itself, if instead of the Schrodinger operator on 
Z d we consider the Schrodinger operator on an arbitrary combinatorial 
or metric graph, whose "global dimension" is greater than 2; this was 
shown in [11]. Here the global dimension is defined as the exponent D, 
such that the heat kernel P(t; x, y), corresponding to the Laplacian on 
the given graph, satisfies the estimate 

(1.3) ||P(i;.,.)|U« =0(r D/2 ), t^oo. 

For the graph Z d one has D = d. 

For the discrete case, in view of (1.2), even a single example of a 
potential such that iV_(— A — aV) x a d l 2 had been unknown for some 
time. Only in [10] the authors suggested a general scheme that al- 
lowed us (for G = Z d ) to construct discrete potentials, such that the 
function N_(— A — aV) has any prescribed asymptotic behavior as 
a — > oo, under very mild regularity conditions. This includes also the 
case iV_(— A — aV) ~ Ca d ^ 2 with C > 0. This scheme was based upon 
some asymptotic estimates for iV_(— A — aV) whose nature differs dras- 
tically from the estimate (1.1) and from the majority of other estimates 
described in [9, 11]: namely, these latter ones are permutation-invariant 
with respect to the potential V. In contrast, the new estimates and 
asymptotic formulas obtained in [10] depend not only on the size of V 
but also on the geometry of its support. More exactly, the potential 
V is supposed to be supported on a sequence of very sparsely placed 
points. The eigenvalue distribution results obtained for such sparse 
potentials are quite flexible. On the other hand, this class of potentials 
is rather special. 

In the present paper we show that the construction of sparse poten- 
tials with a prescribed eigenvalues behavior is not restricted to G — 
Z d , d > 2, but can be, with minor modifications, extended to a wide 
class of combinatorial graphs with D > 2. This we do in Section 2. 

The case D = 2 is critical in this topic, the general scheme of the 
spectral analysis of the Schrodinger operator breaks down at many 
points (as it also happens to the classical Schrodinger operator in M 2 ). 
In Section 3 we show, however, that the approach based upon sparse 
potentials works for the important case G = Z 2 as well (although using 
some more specific analytical tools.) 

Finally, in Section 4 we show that this approach can be successfully 
modified to the Schrodinger operator on such metric graphs, for which 
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the associated weighted combinatorial graph admits the construction 
of sparse potentials. 

Sparse potentials, in the spectral theory of the one-dimensional Schro- 
dinger operator, were introduced by D.Pearson, [8]. Sparse potentials 
in higher dimensions were considered by S.Molchanov and B.Vainberg 
in [6], [7]. Compared with our papers, these authors analyzed prob- 
lems of a different type, and they used another definition of sparseness. 
They did not discuss sparse potentials on graphs different from Z d . 

2. Operators on combinatorial graphs 

2.1. Operator By iG . Let G be a combinatorial graph, with the set of 
vertices V = V(G) and the set of edges £ = E(G). We write e = (v , v') 
for the edge whose endpoints are the vertices v, v'. We then also say 
that v,v' are neighboring vertices, and write v ~ v'. For simplicity, we 
assume that the graph is connected and has no loops, multiple edges, or 
vertices of degree one. We suppose that #V = oo, and that degrees of 
all vertices are finite. On V we consider the standard counting measure 
a. So, a(v) = 1 for each vertex v. Our basic Hilbert space is £ 2 = 
P(G) = L 2 (V; a). We also denote £ q = £ q (G) = L q (V; a), < q < oo. 

With each edge e G £ we associate a weight g e > 0. On £ 2 (G) we 
consider the quadratic form 

(2-1) Bo\f\= E 9e\f(v)-f(v')\ 2 , 

e(z£;e=(v,v') 

with the natural domain 

Doma G = {/Gf 2 (G):a G [/]<oo}. 

This quadratic form is non-negative and closed; recall that by definition 
the latter means that Doma G is complete with respect to the norm 
(a G [/] + H/H 2 ) 1 / 2 . Denote by — A G the self-adjoint operator on £ 2 , 
associated with the quadratic form a G . It is easy to see that on its 
domain the operator A G acts according to the formula 

(2.2) (A G f)(v) = £ f(v')g {v , v , } - f(v) £ g (v , v , } . 

If the weights g e and the degrees degf of all vertices are uniformly 
bounded, then the quadratic form (2.1) and, hence, the operator A G 
are bounded. 

Due to the embeddings l l C I 2 C £°°, any operator, bounded in £ 2 , is 
bounded also as acting from £ l to £°°, with the corresponding estimate 
for the norms. This applies, in particular, to the operators exp(A G t) 
and, therefore, the heat kernel satisfies the estimate \\P(t; ., .)IU°° ^ 
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C < oo, Vt G 1. Our main assumption here is that (1.3) is satisfied 
with some D > 2. Then by the Varopoulos theory, see [13], we have 

(2.3) \\f\\l < Ca G [/], p = p(D) = D(D - 2)"\ V/ G Doma G . 

In particular, (2.3) is valid for any function with finite support. Denote 
by "K{G) the completion of the set of all such functions in the metric 
generated by the quadratic form a G . This is a Hilbert space; the esti- 
mate (2.3) implies that "K{G) can be realized as a space of functions, 
embedded into £ p . It follows from here that any function / G 3i(G) 
tends to zero "at infinity" (notation / — > 0). More exactly, / — > 
means that for any e > the set {v G V : |/(f)| > e} is finite. 

Below we will systematically use the Birman - Schwinger operator 
By iG that corresponds to the Schrodinger operator —A — aV, see, e.g., 
[9] and the references therein. Recall that in our case By jG is the 
operator in the space JC(G), generated by the quadratic form 

(2-4) b v , G [f] = J2v(v)\f(v)\ 2 . 

vev 

The eigenvalue distribution characteristics for the operators —A G —aV 
in £ 2 (G) and By iG in !K(G) are connected by the famous Birman - 
Schwinger principle: for any a > one has 

(2.5) N.(-A G - aV) = nfcT 1 , B VtG ). 

Here n(s, T) = #{n G N : A n (T) > s}, s > 0, stands for the eigenvalue 
distribution function of a compact non-negative operator T. Moreover, 
the next two properties are equivalent: 

(A) The operator By G is compact; 

(B) N-(-Aa - aV) < oo for all a > 0. 

The equality (2.5) allows one to formulate spectral estimates and 
asymptotic formulas for — A G — aV in the terms of the operator By G . 
For instance, the estimate (1.1) means that for V G £ d ^ 2 (Z d ) the oper- 
ator By )G belongs to the "weak Neumann - Schatten ideal" S^/ 2 , and 
that its norm in this ideal satisfies the estimate 

||By, G || Sd/2 < C\\V\\ d/2 ; 

see [1], §11.6, for the definition and basic properties of the classes E g . 

Below we describe the estimates of a completely different nature. 
They are valid for a rather wide class of graphs and for a special, 
but rather restricted class of potentials. To describe these estimates, 
we need some properties of the Green function for the Laplacian A G . 
Recall that we always assume D > 2. 
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We denote by (., .)jf(G) the scalar product in Ji(G), so that 

if,g)x(G) = ac[f,g]- 

Up to the end of this section, we drop the subindex !K(G) in this 
notation and in the notation of the corresponding norm. 

2.2. The Green function. For a fixed v G V consider the linear func- 
tional on !K(G): 

Mf) = f(v). 

In view of (2.3) this functional is continuous on "K(G), and therefore 
there exists a unique function h v G Ji{G) such that 

(2.6) (f,h v ) = f(v), V/eW(G). 
Thus the function h v satisfies the equation 

(2.7) -Ah v (w) = 5 v w 

and is, obviously, real- valued. It is natural to call it the Green function 
for the operator A G . Taking / = h w in (2.6), we obtain 

(2.8) (h w , h v ) = h w (v) = h v (w), Vi>, w G V. 
By (2.3), this implies that 

(2.9) h w (v) — > as w is fixed. 
For w = v the equality (2.8) gives 

(h v , h v ) = h v (v), 

whence all the numbers \j? v = h v (v) are positive. We normalize (in 
"K{G)) the functions h v , i.e., we set h v = fi^hy. 

The function ^(v)" 1 ^ can be also defined as the unique function 
/ G "K{G) that minimizes the quadratic functional ||/|| 2 under the 
condition f(v) = 1. Note also that the number //„ is nothing but the 
capacity of the one-point set {v} with respect to the quadratic form 
a G . 



2.3. Sparse potentials. For a graph G, we introduce the notion of a 
sparse subset Y C V in the same way as this was done for G = Z d in 
[10], Section 6.2. Namely, let Y C G be an infinite set, and let 0~C Y 
stand for the subspace in "K(G) , spanned by the functions {h y : y G 
Y}. We say that Y is sparse if in "K(G) there exists a compact linear 
operator T, such that the operator I — T has bounded inverse and 
the functions e v — (I — T) -1 /^ form an orthonormal (not necessarily 
complete) system in "K{G). We do not discuss here the weakly sparse 
subsets, also introduced in [10]. 
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We cannot claim that every graph G with D > 2 contains a sparse 
subset. In Subsection 2.4 we shall give a simple condition of a geometric 
nature, that guarantees existence of such subsets. 

We say that a function V on G (the discrete potential) is sparse, if 
its support Yy = {v e V : V(v) ^ 0} is a sparse subset. Due to this 
definition, the quadratic form (2.4) can be re-written as 

WAf] = E n«)l(/A)l 2 = E v(v)\(f,(i-T)e v )\ 2 . 

veY v veY v 

If, in addition, V — > 0, we always enumerate the points v EYy, v — v n , 
in such a way that the sequence V(v n ) is monotone. 

Due to the sparseness of Yy, the system of elements h v , v G Yy, is 
close to an orthonormal system in IK(Cr). Hence, the properties of the 
operator By G should be close to those of the diagonal operator with 
the entries V{v), v £ Yy. More precisely, consider the operator 

N v = E Vv^)(-,e v )e v . 

We have 

N V (I-T*)= E VW)i-X)e V) 

veY v 

whence 

l|Nv(I " T*)/||^ (G) = E y^)\U~K)\li {G) = by G [f]- 
This means that 

(2.10) B y , G = (I-T)N^(I-T*). 

The first property of Hv,g follows from (2.10) immediately, cf. The- 
orems 6.3 and 6.4 in [10]. 

Theorem 2.1. Let V > be a sparse potential on a weighted graph 
G with global dimension D > 2. Then the corresponding Birman- 
Schwinger operator By,G is bounded if and only if the function V is 
bounded, and is compact if and only ifV—>0. 
Moreover, 

CH^Iloo < ||B ViG || < C'\\V\U 

where C = ||(I - T) _1 || -2 and C' = ||I - T|| 2 . In the case of compact- 
ness, the following two-sided estimate, with the same constants C, C' , 
is valid for the eigenvalues X n (Byo): 

CV n < \ n (Ky G ) < C'V n . 
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Recall that in the case of compactness we have, just due to the way 
of enumeration, V n \ 0. 

The next result is more advanced. It extends Theorem 6.6 in [10] 
(where it was G = Z d ) to the case of general graphs. The proof survives. 
Still, we reproduce it here, to make it possible to read the present paper 
independently. 

We derive this result under an additional condition on the potential 
V: we assume that the sequence {V n } is moderately varying. This 
means that V n \ and V n+ i/V n — > 1. We shall use a result of M.G. 
Krein, see Theorem 5.11.3 in [3]. Below we present its formulation, 
restricting ourselves to the situation we need. 

Proposition 2.2. Let H > and S be self-adjoint and compact oper- 
ators. Suppose rankH = oo and the sequence of non-zero eigenvalues 
A n (H) is moderately varying. Then for the operator M = H(I + S)H 
one has 

A n (M) ~ A£(H). 

Now we are in a position to prove the following result. 

Theorem 2.3. Let V > be a sparse potential, such that the numbers 
V n form a moderately varying sequence. Then 

A n (By G ) ~ V n . 

Proof. In view of (2.10) the non-zero spectrum of the operator By G 
coincides with that of the operator 

(2.11) My := N y (I - T*)(I - T)Ny = Ny(I + S)Ny, 

where 

Now we apply Proposition 2.2 to the operators H = Ny and My, 
given by (2.2) and (2.11) respectively. All the assumptions of Proposi- 
tion are evidently satisfied, and we get the desired result. □ 

In particular, taking V(v n ) = n~ 2 / D , we obtain a potential V such 
that the eigenvalues A n (Bycr) asymptotically behave as n~ 2 l D '. This 
solves a problem discussed in Section 1. 

2.4. On the existence of sparse subsets. Given a vertex v E V, 
let us consider the function /„ that vanishes at all vertices w ^ v and 
f v (v) = 1. By the extremal property of the Green function h v , we have 



= iOlM 2 < Wfvf = ^a[fv] = ^9{v. 
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whence 



l4 = K{v) > ^2g( v ,w) 




Given a number R > 0, we say that a vertex v R-mild if 



9{v,w) < -R- 



Proposition 2.4. Suppose that a graph G is such that D > 2 and £/ia£ ; 
for some R > 0, G contains infinitely many R-mild vertices. Then G 
contains a sparse subset. 

Proof. We fix an infinite subset G' C G, consisting of i?-mild vertices. 
Choose any double sequence {e mn }, m,n G N; m ^ n, of positive 
numbers, such that e mn = e nm and J2 £ mn < R~ 2 - Suppose that the 
vertices v±, . . . ,v n -i are already chosen (here v\ is arbitrary). Then, 
in view of (2.9), we can choose a point v n G G' in such a way that 
\(h Vk ,h Vn )\ < Ekn for k — 1, . . . , n — 1. As a result of this induc- 
tive procedure, we get a sequence of vertices {v n } C G', such that 
( h v n ,h Vn ) = /4 n and l(^m>^„)l < £ mn if m < n. By (2.8), the same 
inequality holds for m > n. Since all the vertices v n are i?-mild, we 
have \(h Vn ,h Vn )\ < Re mn . 

It follows that the Gram matrix M = {(h Vn , h Vn )} is such that M — I 
is Hilbert-Schmidt and, moreover, ||M — 1||© 2 < 1. By Theorem VI. 3. 3 
in [3], this implies the sparseness of the set Y, with T G ©2; see 
[10], section 6.2 (where the case of G = Z d was analyzed) for more 
details. □ 



It is well-known that a lot of complications arise in the study of the 
classical Schrodinger operator in R 2 , compared to the higher-dimensio- 
nal case. They are related to the non-signdefiniteness of the funda- 
mental solution, its logarithmic singularity, the absence of the proper 
sharp Sobolev embedding theorem, and so on. One of the important 
consequences here is the absence of a direct analogy of the Birman- 
Schwinger principle (2.5), since the closure of the space of compactly 
supported functions in the metric of the Dirichlet integral is not a space 
of functions any more. 

We encounter the same complication if we try to extend the general 
reasoning in the previous section to the case of operators on combina- 
torial graphs with global dimension D = 2. Similar to the case of 1R 2 , 
much more specific instruments are needed. 



3. Operators on Z 2 
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In this section we present an approach to the eigenvalue analysis of 
the Schrodinger operator on Z 2 . 

We consider the graph G whose set of vertices V is the lattice Z 2 
which we understand as embedded in the natural way into M 2 . Two 
vertices are said to be connected with an edge if they are the closest 
neighbors in R 2 , and the weight of any such edge is taken as 1. It is 
convenient to identify a vertex v G Z 2 with its Cartesian co-ordinates 
x = ( Xl ,x 2 ) e Z 2 C R 2 . 

3.1. The Birman-Schwinger principle, a Hardy type inequality, 
and the space J£ . Following Section 2, we introduce the quadratic 
form 

(3-1) M/]= £ \f(x)-f(x')\ 2 . 

x,x'£Z 2 ,x~x' 

This quadratic form is bounded in I 2 = I 2 (7?) and defines the self- 
adjoint operator —A. For a bounded real- valued function V(x) on Z 2 
and a coupling constant a > 0, we consider the Schrodinger operator 
—A — aV, and our concern is in the study of the behavior of i\L (—A — 
aV) as a — > oo. 

So far so good, however if we try to apply the Birman-Schwinger 
principle, as we did it above, we fail. The reason for this lies in the 
fact that the space "K, the closure of the space of finitely supported 
functions in the metric (3.1), is not a space of functions any more. In 
fact, it is easy to construct a family of functions u n on Z 2 , converging 
to a nonzero constant point-wise, so that a^f-u^] — > 0. This is also 
related to the behavior of the heat kernel for —A on Z 2 . Standard 
calculations by means of the Fourier series show that the heat kernel 
P(t; ., .) decays exactly as t^ 1 as t — )■ oo, so the global dimension equals 
2 and the Varopoulos theory does not apply. 

In order to handle these inconveniences, it is sufficient to impose the 
Dirichlet boundary condition at one point, say at x — = (0, 0). 

So, denote by i\ the (obviously closed) subspace in £ 2 consisting 
of functions with zero value at 0. The quadratic form (3.1) defines 
a self-adjoint operator — A in £q, and we denote by — A — aV the 
corresponding Schrodinger operator. 

Since the subspace l\ has codimension 1 in I 2 , the standard applica- 
tion of the variational principle gives 

(3.2) iV_(-A - aV) < AL(-A - aV) < AL(-A - aV) + 1, 
therefore, for our aims, it suffices to study the operator — A — aV. 
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The minor reduction above brings a great advantage: applying the 
Birman-Schwinger principle becomes possible now. To explain this, we 
start by establishing a Hardy type inequality in l\. 

Proposition 3.1. For some constant C and for any function f e t\ 
with compact support, 

(3.3) \f^)\ 2 \x\- 2 (\og\x\+2)- 2 <Ca z2 [f]. 

xei?\{o) 

Proof. The inequality (3.3) is proved in a way, similar to the discrete 
Hardy type inequality in Sect. 4 in [10], by interpolating a function on Z 2 
to a function on IR 2 . In all squares in M 2 with vertices in Z 2 , we apply 
the same piece- wise linear interpolation as in [10]. Additionally, for 
the four central squares, the ones having the vertex 0, we multiply the 
interpolating function by the cut-off 9?(|a;|), which equals for \x\ < \ 
and equals 1 for \x\ > |. The inequality (3.3) follows now by applying 
the classical Hardy inequality in R 2 with log term (valid for functions 
vanishing in a fixed neighborhood of the origin) to the interpolant. □ 

We define the space "K as the closure of the set of finitely supported 
functions in i\ in the metric (3.1). It follows from (3.3), that on this 
set the convergence with respect to (3.1) implies the convergence in £ 2 
with the weight as in (3.3), therefore is a space of functions. 

Now we are able to define the Birman-Schwinger operator in 
the space CK by means of the quadratic form (2.4), obviously, closed. 
By the general Birman-Schwinger principle and (3.2), we have 

n(a-\ B^ 2 ) < AL(-A - aV) < n(a~\ B^ 2 ) + 1. 

3.2. The Green function. It follows from (3.3) that the linear func- 
tional (fi v (f) = f(x) is continuous in 0~C and defines a unique function 
h x G "Ko such that 

(3.4) (f,h x ) = f(x), V/GMfl. 

Thus the function h x satisfies the equation (2.7); it is real- valued and 
will be called, similar to the considerations in Section 2, the Green 
function for the operator A on Z 2 . The relation (2.8) holds too, but 
at this point the analogy with Section 2 stops. The reason for this 
is that we do not have the inequality (2.3) any more, and the Hardy 
inequality (3.3) does not imply that h x (y) — > as \x — y\ — > oo. More- 
over, the asymptotic formula (3.9) below shows that in Z 2 -setting this 
property fails. Nevertheless, we define, again, h x = n x l h Xl where 
/i 2 = (h x ,h x ) = h x (x). Again, as it was in the case of graphs with 
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D > 2 (see the end of Subsection 2.2), the function n~ l h x is the unique 
minimizer (in the space < Kq('L 2 )) of the functional (3.1) under the con- 
dition f{x) = 1. 

We are going to find an explicit expression for h x . To do this, we 
consider first the fundamental solution for the Laplacian on Z 2 without 
the boundary condition at 0, i.e., the function Q(x), x G Z 2 , such that 
AS(x) = 5q. Such function, if exists, is not unique, it is defined up 
to Z 2 -harmonic additive term. We are going to deal with a specific 
function 9- 

For Z d with d > 2 such solution can be easily found by means of the 
Fourier series, see the formula (6.1) in [10]. In the 2-dimensional case, 
however, the integral in [10], (6.1), diverges, so a different approach is 
needed. 

The first formula for 3(x) was found in [4]. We do not reproduce it 
since it is not convenient for further calculations. We are going to use 
another formula for the fundamental solution, presented, for example, 
in the book [12], formula (1) in Chapter 15: 

(3.5) 9<*) = <ao-» / ^^de, 

(— 7T,7T) 2 

where 

Z(9) = Z{9 U 9 2 ) = 1 - -(cos 0i + cos# 2 ), x9 = x 1 9 1 + x 2 9 2 . 
For our aims it is more convenient to use the exponential form of (3.5): 

(3.6) 9(*) = (*)-" / ^f|^». 

( — TV, 1l) 2 

So, 9(0) = 0; it is also known, see [12], Chapter 15, that S(x) has 
the asymptotics 

(3.7) 5(x) ~ -(21og|x| +log8 + 2 7 ), | x\ — > oo, 

uniformly in all directions (7 is the Euler constant). 
Now we construct the function H x (y), as follows: 

(3.8) H x (y) = -9(x -y) + 9(x) + 9(y). 

One can easily check that H x (y) = H y (x), H x (0) = and —A y H x (y) = 
b~ y x . It will be shown below that H x = h x . 
By the asymptotics (3.7), 

(3.9) H x (y) ~ i(21og|x|+lo g 8 + 2 7 ), \y\ -> 00 
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2 

(3.10) H x (x) = 2Q(x) ~ — (21og |a?| + log8 + 2 7 ), \x\ ->■ oo 

The following property is the most important one. 
Proposition 3.2. For any x E Z 2 \ 0, the function H x (y) belongs to 

Proof. We denote 1 2 = (0, 1), fix x E Z 2 \0 and consider the expression 
tfeifsfe) = H x (y + 1 2 ) - if x (y). By (3.6), 

d 2 H x (y) = (2tt)- 2 J ( e <(»+ia)« - e *« - e i(i/-*+ia)« + e %-*)0) ^ 

( — 7T,7r) 2 

* 1 eMwO)de. 

(-7T,7r) 2 

Considered as a function on the torus (— 7r,7r) 2 , the denominator in 
(3.2) has the only zero at the point 9 = 0, and Z{6) x \9\ 2 near = 0. 
Both functions e a2d — 1, e - " 61 — 1 vanish at 9 = 0, the main parts being 
degree one homogeneous in 9. Therefore the function 

( e « a g-l)(e-fa*-l) 

* (£ ° = Z(0) 

is bounded on (— n,n) 2 , in particular, it belongs to L 2 ((— ir, tt) 2 ). So, 
the expression (3.2) represents the Fourier coefficient of the function 
K(9), and by the Plancherel equality, we have 

d 2 H x (-)E£ 2 (Z 2 ). 

In a similar way, the function d\H x {y) = H x (y + li) — H x (y) (where 
li = (1, 0)) belongs to i 2 (Z 2 ). These two facts mean that a Z 2[H x ] < oo 
for any x. 

It remains to show that the function H x can be approximated in the 
metric a G by finitely supported functions. This is done in the standard 
way, by setting H x m) (y) = H x (y)F m (\y\), where F m (s) = l,s < 1, 
F m (s) = (logm) _1 (logm — logs),s E (l,m), and F m (s) = 0, s > m. 
Using that H x is bounded and the above properties of djH x , it is easy 
to show that a^2 [H x — H x "^] — > as m — > oo. □ 

So, the functions H x and h x , both in "Kq, are the solutions of the 
same equation. Due to the uniqueness in the construction of h x , we 
have 

hx H x , 
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and all the properties established in this section for the functions H x 
hold now for the functions h x . 

3.3. Sparse sets in Z 2 . The definition of a sparse set given in Sub- 
section 2.3, requires a minor modification in the case of Z 2 : the space 
!K should be replaced by JC - The next statement is an analogue of 
Proposition 2.4. However, the proof is a little bit different, it uses the 
asymptotic formula (3.7) and its consequences. 

Proposition 3.3. The graph Z 2 contains a sparse set. 

Proof. For the functions h x introduced in Subsection 3.2, we have 

(h x ,h y ) = h x (x)~2h y (y)-2h x (y). 

For a fixed x G Z 2 , we use the asymptotics (3.9) and (3.10) as y — > oo. 
This gives 

(3.12) (h x J y ) ~ ^(log |y|)-2(21og \x\ + log8 + 2 7 ). 

It is clear from (3.12) that (h x , h y ) — » as x is fixed and y — > oo. 

The rest of the proof uses the same inductive procedure as in Subsec- 
tion 2.4. Fix a double sequence {e mn }, m,n G N; m ^ n, of positive 
numbers, such that e mn = e nm and £ mn < 1- Suppose that the points 
x±, . . . ,x n _i are already chosen {x\ is arbitrary). In view of (3.12), we 
can choose a point x n G Z 2 \ such that 

(3.13) \(h Xm ,h Xn )\ < e mn , Vm < n. 

As a result of this inductive procedure, we obtain an infinite sequence 
x n G Z 2 \ such that (3.13) is valid for all m,n G N; m < n. By 
symmetry, this inequality is satisfied also for m > n. 

It follows (exactly as it was in Subsection 2.4) that the Gram matrix 
M = {{h Xm) h Xn )} is such that ||M — I||e 2 < 1, and hence, by Theorem 
VI. 3. 3 in [3], the set Y = {x n } is sparse. □ 

Now we define a sparse potential on Z 2 as such whose support is a 
sparse subset. The spectral properties of the corresponding operator 
and hence, those of the operator — A — aV, are studied in the 
same way as this was done in Section 2. This leads to the analogues of 
Theorems 2.1 and 2.3. The next statement is a combination of these 
both results. The proof remains the same and we skip it. We only 
remind that the operator T in the formulation is the one, appearing in 
the definition of a sparse subset, cf. Subsection 2.3. 

Theorem 3.4. Let V > be a sparse potential on Z 2 . Then 
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1. The corresponding Birman-Schwinger operator ~B> VI 2 is bounded 
if and only if the function V is bounded, and is compact if and only if 
V ->■ 0. 

Moreover, 

C||V^||^oo < 1 1 By J2 || 5: C ||V||#x>, 

where C = ||(I — T) 1 1 1 — 2 and C = ||I — T|| 2 . In the case of compact- 
ness, the following two-sided estimate, with the same constants C, C , 
is valid for the eigenvalues A„(By^ 2 ): 

CV n < A n (B^ 2 ) < C'V n . 

2. If, in addition, the numbers V n form a moderately varying se- 
quence, then 

^n(By iZ2 ) ~ V n . 
4. OPERATORS ON METRIC GRAPHS 

Let T be a metric graph, with edge lengths l e . As in [11], we associate 
with T the combinatorial graph G = G(T), with the same sets V, £ and 
the same connection relations, and we suppose the general conditions 
on G formulated in the beginning of Section 2 to be fulfilled. To any 
edge e of G(T) we assign the weight g e = I" 1 . The basic Hilbert space 
is now L 2 (T), with respect to the measure induced by the Lebesgue 
measure on the edges. 

The Sobolev space JC 1 (r) consists of all continuous functions ip on 
T, such that ip G H l (e) on each edge and J r (|<^'| 2 + |y?| 2 )rf:r < oo. The 
operator —A on Y is defined via its quadratic form ar[</?] = f r (\ip'\ 2 dx, 
with the form-domain ^(T). The Laplacian A acts as y?" on each 
edge, and the functions from DomA satisfy the natural (Kirchhoff) 
conditions at each vertex. 

Here we list our assumptions about the graph T. 

Assumption 4.1. One of the following holds: 

(i) The edge lengths l e of the graph V are bounded from above, the 
corresponding combinatorial graph G(T) has global dimension D > 2 
and contains a sparse subset, say Y . 

(ii) The graph G(r) is Z 2 . 

To define the Birman-Schwinger operator, we need to consider two 
cases of Assumption 4.1 separately. 

In the case (i), by theorem 4.1 in [11], the graph V has the same 
global dimension D. Therefore the inequality of the form (2.3), with 
a r replacing a G , holds for all functions in ^(T) with compact support, 
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and consequently the space IK(r), the closure of the set of compactly 
supported functions in the metric a r , is a space of functions (see [11], 
where the details are presented.) 

In the case (ii), i.e., G(T) = Z 2 , one should consider the subspace 
!Kg(r) in [K 1 (r), consisting of functions, vanishing at the vertex 0. A 
Hardy type inequality, derived easily from (3.3), implies that the space 
J-Cq(T), the closure of the set of compactly supported functions in ^Kq(T) 
in the metric ap, is a space of functions as well. 

Further on, to unify notations, we suppress the subscript when 
dealing with JC (r), so, for G = Z 2 , JC(r) means IKo(r). 

An analogue of the quadratic form (2.4) is given by 

(4.1) b v M = jv\y\ 2 dx, 

where V > is a function in L 1 (r). We denote by By,r the operator 
in "K(T), generated by the quadratic form (4.1). 

We will construct potentials V on T, whose support is a vicinity (in 
T) of the set Y . By using such potentials, we will be able to prove the 
following statements that can be considered as analogues of theorems 
2.1 and 2.3. 

Theorem 4.2. Let for a graph Y Assumption 4-1 be satisfied. Let a 
sequence p n \ , such that n}l 2 p n — > 0, be given. Then there exists a 
function V G L 1 ^), V > 0, such that 

C Pn < A n (B y , r ) < C' Pn , Vn G N, 

where the constants C , C G (0, oo) do not depend on the sequence {p n }. 

Theorem 4.3. Suppose in addition that the sequence p n is moderately 
varying. Then the function V G -^ 1 (r), V > 0, can be chosen in such 
a way that 

X n (B V:T ) ~ Pn . 

Below we restrict ourselves to the proof of Theorem 4.3. The proof 
of Theorem 4.2 is much simpler and can be easily reconstructed by the 
analogy. As in [11], section 4.3, for the analysis of the operators By r 
we use of the orthogonal (in Ji(T)) decomposition 

(4.2) ft(r) = JC p i © -k v . 

Here the subspace "K p \ is formed by functions tp G !K(r), which are 
linear on each edge e, and IK© is formed by functions vanishing at each 
vertex v. Below we denote by By^pi, By,r ; D the operators in these 
subspaces, generated by the quadratic forms byr [ ^ P i and byr [ "Kv 
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respectively. Given a function ip G Jt(T), we denote by ip p \, ipv its 
components in the decomposition (4.2). 

Proof of Theorem 4-3. Let V be a sparse potential on the combinatorial 
graph G(T), supported on a sparse set Y G V, such that V(v n ) — p n \ 
0, i> n G V. Along with the operator B y G( - r ^ let us consider the operator 

By r in the space IK(r), associated with the quadratic form 
This quadratic form vanishes on Ji®, and 

(4.3) by^]=hy MT) [ V \V]. 

It follows that the non-zero spectra of the operators B y r and B y G ^ 
coincide. So, if we allow such point-supported potentials, then all the 
results of Section 2.3 carry over to the metric graphs. 

Now we show how to construct a "genuine" potential V G L 1 ^) that 
possesses the desired properties. Suppose for simplicity that Y = Y y 
contains no neighboring vertices. This can be always achieved by a 
thinning down the original sparse set. Choose a sequence e n > 0, such 
that e n < min/ e . Further assumptions about the behavior of e n will be 

eBv n 

imposed later. For the calculations below, it is convenient to denote 
by £ the set of all edges, such that one of its ends lies in Y. 
Define now the potential F on T, as follows: 



(4.4) V{y) 



£n ife^ if dist(y,v n ) < e n for some n, 
otherwise. 



We have 

(4.5) by,r[<p] = by, r [</?pi] + by,r[<pi>] + 2 Rebv,r[Vpi, <Pv\- 
We shall inspect each term separately. 

Any edge e G £ can be written as e = (v n ,v' n ) where v n G Y and 
v' n Y. The vertex v' n that corresponds to the vertex v n is not unique, 
so that here we have some slight abuse of notation, however what is 
important is that v' n is determined uniquely by e. We identify any such 
edge with the interval (0, l e ), so for y G e we have 

VAV) = C 1 (f( V n)(le ~y) + fW n )y) ■ 

From here we derive 

by,rM = Yl ~nTzf / \f( v n)( l e ~y) + ¥K)y\ 2 dy. 

_ t e c n aeg v n Jq 

ee£ 
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In each summand the term that will be shown to be dominant is 

Pn r Sn Pn 

e n deg v n J deg v n 

So, by (4.3) 

= b v,G(r)M- 

To estimate the difference byr^pi] — b^> g/pJ^pi] we use the elemen- 
tary inequality \\A\ 2 - \B\ 2 \ < \A - B\(\A\ + \B\), to obtain 

\\<p(v n )(l e - y) + <p(v' n )y\ 2 - \<p(v n )\ 2 l 2 e \ 

< y\cp(v n ) - tp(v> n )\ ((2l e - y)\<p(v n )\ + y\<p(i/ n )\) < 2l e y{\ V {v n )\ 2 + |y?«)| 2 ). 

Now the summation over all edges e G £ gives 

(4.6) |by,rb P i] - by iG(r) [</? p i]| < Yl R{e,e n ;<^\V) 

e=(v n ,v'„)€£ 

where 

R(e,e n ) = -^^Mv n )\ 2 + Hv' a )\ 2 ). 

Choosing e n decaying fast enough, we can grant an arbitrarily fast 
decay of the eigenvalues of the operator corresponding to the quadratic 
form in the right-hand side of (4.6). This yields that 

An(By i r ;p i) ~ p n . 

Now we go over to the quadratic form by[</?i)]. The associated opera- 
tor By r;!) can be identified with the direct orthogonal sum of operators 
By e , e G £ acting in the spaces if 1,0 (e), and the Rayleigh quotient for 
each operator By e is 

' i.y(vMv)\^ ueHl ,« (e) 

where the weight function V(y) is defined by (4.4). For the eigenvalue 
estimates of such operators we use the following 

Proposition 4.4. Let V{y) > be a monotone function on an interval 
e = (0, a), a < oo, such that j e ^JV{y)dy < oo. Then for any A > 



n(A,By e ) < -A" 1 / 2 / ^W)dy- 

K Je 



Proof. For a = oo, due to the Birman - Schwinger principle, this is an 
equivalent reformulation of the well known Calogero estimate ([2]; see 
also [5], Theorem XIII. 9(b)). The case a < oo follows from here by the 
standard variational argument. □ 
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Applying Proposition 4.4 to each edge e = (v n ,v' n ) G £, we obtain 
(4.7) n(A,B v , ra ) = < ^" 1/2 E (j^ V2 



ee£ e<s£ 



If the series on the right converges, then n(A, B^d) = 0(A -1 / 2 ) (or, 
equivalently, \j(Bv,r;v) < Cj~ 2 ), and hence, the operator By ; r ; D does 
not contribute to the asymptotics of any greater order. 

It remains for us to inspect the last term in (4.5). It is equal to 

%>]=2Re J2 WT^T [ £ \v(vn)(le-y)+rtv' n )y)M^dy. 

_ t e c n aeg v n Jq 

e=(vn,v'„)eE 

The main contribution to this sum is given by the similar expression, 
say fiolv 9 ]) with the first factor in the integrand replaced by l e tp(v n ). 
For the spectral estimates we choose one more sequence j n , vanishing 
as n — > oo. 

Each term in f2oM does not exceed 

''" \<PM\ [ \(f V (y)\dy 



l e e n deg v n 



< lnPn\v( V n)\ 2 + -. 7T~ ^ £ nl l^Mf ''dy ■ 

A-f n (degv n y J 

Here the first terms on the right correspond to the quadratic form that 
by Theorem 2.1 generates an operator whose eigenvalues, due to the 
factor 7„, behave as o{p n ). The second terms are quadratic forms of 
the same type as the ones considered above, for the operator By 5 r ; D, 
with the same weight function. Applying the estimate (4.7) (and again, 
choosing e n decaying fast enough), we see that the corresponding esti- 
mate is of the Weyl type, and so, this terms do not affect the spectral 
asymptotics of any order greater than 0(A -1 / 2 ). 

The estimates for the quadratic form £l[ip] — ^oM are much easier. 
We have 

|nM-noM|<2 £ PnWK)-^v n )\ r, ( )|d 

^ Ij deg v n J 

Each term here can be estimated through 

ti(\<PM)\ 2 + \<p(vn)\ 2 )e n 



+ / \<Pv{y)\ 2 dy. 
Jo 



The first terms here give a quadratic form whose eigenvalues decay as 
fast as we wish, provided the sequence e n is chosen in an appropriate 
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way. The same is true for the second terms, after we apply the estimate 
(4.7). □ 
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